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Работа посвящена описанию модифицированного метода нечеткой кластеризации полутоновых 
изображений, который на каждой итерации выполняет динамическое преобразование исходных данных на 
основе сингулярного разложения с автоматическим отбором наиболее значимых столбцов матрицы левых 
сингулярных векторов. Это дает возможность повышения чувствительности сегментации. Предложенный метод 
показан на примере алгоритма нейро-фаззи кластеризации $ЕСМ. Результаты экспериментальных исследований 
были получены при обработке реальных полутоновых медицинских изображений. При этом было достигнуто 
лучшее выделение объектов интереса и структуры изображения в целом. 
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Введение 

Существует множество различных 
алгоритмов кластеризации, основанных на 
детерминированных, статистических, не- 
четких и эвристических моделях. Фунда- 
ментальное различие между типами этих 
алгоритмов существенно зависит от харак- 
теристик доступных данных. 

Кластеризация изображений является 
сложной процедурой, способ выполнения 
которой существенно зависит как от их 
характеристик, так и от особенностей реша- 
емой задачи [1]. Разнообразие целей обра- 
ботки и типов изображений, различные 
недостатки существующих методов (напри- 
мер, недостаточная чувствительность при 
обработке низкоконтрастных изображений, 
большое количество управляющих парамет- 
ров, низкая производительность и т.п.), 
зависимость эффективности метода класте- 
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ризации от исходных данных усложняет 
выбор подходящего алгоритма при решении 
конкретной задачи. 

Таким образом, разработка новых/ 
модификация существующих методов 
кластеризации остается актуальной задачей. 

Постановка задачи 

В настоящее время для решения 
задачи кластеризации часто применяются 
нечеткие алгоритмы, нейросети, а также 
нейро-фаззи методы, позволяющие объе- 
динить достоинства обоих подходов. В 
работах [2, 3, 4] был описан алгоритм 
$ЕСМ, представляющий собой объедине- 
ние сети Кохонена и алгоритма нечеткой 
кластеризации ЕСМ. При обработке 
изображений с помощью алгоритма $ЕСМ 
существует возможность выполнения 
динамического преобразования функции 
принадлежности (на основе применения 
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методов повышения контраста), что позво- 
ляет повысить чувствительность кластери- 
зации. Выбор метода динамической 
обработки оказывает существенное влия- 
ние на конечный результат. 

В работе [5] было предложено 
применение метода сингулярного разложе- 
ния при визуализации результатов нечет- 
кой кластеризации с использованием 
композиции исходного изображения и 
нечеткой функции принадлежности, что 
позволяет обеспечить повышение детали- 
зации в областях интереса при выполнении 
сегментации. 

В данной работе предложен модифи- 
цированный алгоритм $ЕСМ для кластери- 
зации полутоновых изображений, который, 
за счет выполнения динамического преоб- 
разования с применением сингулярного 
разложения исходных данных на каждой 
итерации, позволяет обеспечить повыше- 
ние чувствительности сегментации. 

Изложение основного материала 

Предложенный в данной работе 
алгоритм состоит из следующих шагов: 

1. Инициализация начального числа 
кластеров с и значений матрицы цент- 


ров нечетких кластеров У0; экспонен- 
циального веса нечеткой кластеризации 
т; коэффициента увеличения числа 
нечетких кластеров №. 


2. Формирование начальных векторов ве- 
сов нейронов ЗОМ, на основании мат- 
рицы центров нечетких кластеров пре- 


дыдущей итерации Ут ‚ количество 


которых [№, -с,1] (одномерная структу- 
ра — столбец). Дополнительные значе- 


НИЯ Ут вычисляются методом про- 
порционального распределения [3, 4]. 

3. Кластеризация исходных данных с 
помощью ЗОМ, происходящая в два 
этапа: грубая и тонкая настройка весов 
нейронов. При этом выполняется 
динамическое уменьшение количества 


нейронов, описанное в работе [3, 4]. 
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4. Формирование новых значений цент- 


ГА 


ров нечетких кластеров У от 


путем вы- 


бора значимых нейронов, полученных в 
результате обучения ОМ. 
5. Вычисление текущих значений функ- 


ции принадлежности и! 
—2 


с | Б;у т—1 (УК Е Вов 
рые" [и © 
р [= В; УТ Е 1.....п} 
причем матрица 2) 
следующим образом: 


аи) абы). о 


(9 = 1... п], УК = 1...) 


вычисляется 


где матрица А = Г. 

6. Динамическое преобразование исход- 
ных данных на основе сингулярного 
разложения, которое выполняется 
следующим образом: 

6.1 Формируется матрица 1 а размер- 
ностью пх (@ +5), которая содержит 


исходное изображение и матрицу 0. 
Причем п - количество пикселей 
исходного изображения а @ - 
количество цветовых компонент (для 
полутоновых изображений & =1). 


6.2 Полученная матрица р подвергается 
сингулярному разложению [6], в 
результате чего формируются три 


матрицы: / ь (содержит левые сингуляр- 

ные вектора), 5; (содержит собствен- 

ные числа) и У; (содержит правые 
сингулярные вектора) такие, что: 

П=П.5,. (У). 3) 

При этом размерности матриц /[ Ти 

1 ь совпадают, а размерность матриц 5; 

и У, имеет значение (4+ с)х(@+с) 


(используется экономное представле- 
ние сингулярного разложения). 


6.3 Все столбцы матрицы 1 - масштабиру- 
ются на отрезок [01] и интерпрети- 
руются как матрица размерности 
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аухахх(@+с), где 4у и 4 - раз- 
мерность исходного снимка по верти- 
кали и горизонтали, соответственно. 


6.4 На основе матрицы У, вычисляется 


вектор коэффициентов | по 
следующей формуле [5]: 
Ш -+сС + 


С; = У, ). + У): 12 
1 > 5); ра 5) .(4) 
ТЕ [, 4+ с] 
6.5 Вектор коэффициентов С 


упорядочивается по убыванию, а его 
элементы нормируется так, чтобы их 
сумма была равна 1. 

6.6 Формируется вектор АС ‚ содержащий 
разности для каждой пары соседних 
элементов вектора С. 

6.7 Вычисляется значение @С,„ 
формуле: 

+с-1 
>хас; 
ее ‚ (Ч Ста +АСтих) |5) 
2(4+с-1) 4 
где АСии и АЧСах -— мини- 


по 


мальный и максимальный элемен- 

ты вектора 4С ‚ соответственно. 
6.8 На основе значения 4С„ при просмотре 
элементов вектора @С с конца 
выбирается такой индекс Пиах, Для 


которого выполняется условие: 


ас, > аСа . (6) 
6.9 Если выполняется условие: 
И С (7) 


то динамическое — преобразование 
исходных данных не будет более 
осуществляться, и происходит переход 
к пункту 7. Необходимость этого шага 
обусловлена ухудшением детализации 
исходного изображения при получении 
большого числа информативных столб- 


цов матрицы п. Условие (7) было 


получено в ходе проведения 
экспериментов. 
6.10 Элементы вектора АС с индек- 


сами от 1 до ах +1 нормируется так, 
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чтобы их сумма была равна 1. 


6.11 Формируется изображение / ? по 
формуле: 
2 тах +1 ] 
Ту,х ы >. Гух,] ‘АС у, (8) 
1=1 
уЕ [, ау х Е [,ах] 

6.12 Изображение [ 7 масштабируются 
на — отрезок [ОП, после чего 
интерпретируется, как исходные 
данные Х. 

6.13 Выполняются шаги 2 — 5 для 


ГА 
получения обновленных матриц Ухот И 


ГА 
и. 


7. Формируется матрица У т по 


формуле: 
и 7 я И (9) 
(УК {,...с}. М) е ,...4}) 


Эта матрица будет использована в 
начале следующей итерации. 
й 
8. Вычисляются значения А, как среднее 
по матрице расстояний между цент- 


рами У т и Ут Та также критерии 


и, иу' По формулам 


К=И=1 = 
УЕ, == . (10) 
и} 
с п яд 
254.) 
1 _ К=ЦЕ 
Е (р 
где 4 - количество информативных 


признаков каждого вектора исходных 
данных (при сегментации полутоновых 


х е 
изображений равно 1), ии - мини- 


мальное Евклидово расстояние между 
центрами нечетких кластеров. Умень- 
шение показателями Ксие-Биени и 
увеличение показателя нечеткости 
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характеризуют улучшение качества 
нечеткого разбиения на группы. 
9. Если выполняется условие: 


1 тах 
СЕ2Сь ; 


р 
причем значение С» при поиске 


(12) 


конечного результата вычисляется по 
формуле: 


о. 
о А (13) 
где Св” максимальный из коэф- 


фициентов С р ‚ полученных в процес- 


се обучения, то запоминаются следую- 


. лмах _ д шах _ 
щие значения: ЛД,“ =Д,, СЕ =Сь, 
шах _ тах _ 1 
ИЕ, И, ЗУ 


10. Если не выполняется условие: 


1 
Ду, <= или 


(У -У | <виуь УЕ! |< =), (14) 


#1 #1 
где У И У. — показатели Ксие- 


Биени и нечеткости предыдущей 
итераций, соответственно, то переход к 


пункту 2. 
11. Если выполняется условие 


С <Сь^ и (4, > АуйХ или 
(Му <Аухир, >), — (5) 
причем коэффициенты р „И РГ 


определяются выражениями: 


1 шах 
скс“ 1 
Ре 1 тах. -тах ши ’ (16) 
тах(С СЕ ) р -СЕ 
1 
А — АХ | и 
РА, = г тах). | оО] 
тах(А” , ДтаХ ) (тах | _ до 
где С т И до — минимальные 


значения параметров С» и м соот- 


шах 
ветственно, а ( — максимальное 
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значение критерия Д,, то происходит 


возврат к сохраненным значениям 
матриц нечеткой функции принадлеж- 


ности и“ и центров нечетких 
кластеров Ут › которые и являются 


результатом обучения. 

Экспериментальные результаты. 
Представленный алгоритм был применен 
при обработке различных полутоновых 
изображений, примером которых являются 
снимки, представленные на рис. 1, 2. 

На рис. 1 а приведена томограмма 
головного мозга, выполненная с целью 
диагностики наличия гематомы, а также 
определения области ее влияния. На основе 
анализа гистограммы исходного изображе- 
ния (рис. 1 6) этот снимок нельзя отнести к 
классу низкоконтрастных. Однако диагнос- 
тирование гематомы по исходному изобра- 
жению затруднено из-за ее расположения 
на низкоконтрастном участке снимка 
(фрагмент, обведенный прямоугольником), 
а область ее влияния вообще визуально 
неразличима. Введение рентгеноконтраст- 
ного вещества (рис. 1 в) хоть и позволяет 
более четко выделить гематому, однако не 
выявляет области ее влияния. 

На рис. 2 а представлено полутоновое 
медицинское изображение — спин-спиновая 
Т2 релаксация ЯМР участка головного 
мозга. По гистограмме исходного снимка 
(рис. 2 6) можно заключить, что это темное 
низкоконтрастное изображение. 

При проведении экспериментов в 
базовом алгоритме ЕСМ динамическое 
сжатие и обработка нечеткой функции 
принадлежности не использовались. Как 
для базового, так и для модифициро- 
ванного алгоритмов перед применением 
карты Кохонена число нечетких кластеров 
увеличивалось в 4 раза, а также осущест- 
влялось упорядочение центроидов на осно- 
ве матрицы Евклидовых расстояний. Вы- 
бор исходного числа нейронов после 
обучения карты Кохонена осуществлялся 
на основе максимального соответствия. 
Количество нечетких кластеров с=6б в 
обоих случаях. Визуализация результатов 
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нечеткой кластеризации выполнялась на 
основе сравнения с исходными данными [4] 


(Г), рлхе[ 
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08 093 


Рис. 1. Рентгеновская томограмма головного мозга: а — исходный полутоновой снимок 
(204х201); б - его гистограмма; в — результаты введения рентгеноконтрастного вещества 


ИГ), рахе[5 


1200 | 


1000 | 


800 


00 


400 


а 
Рис. 2. Спин-спиновая Т2 релаксация ЯМР участка головного мозга: 
а — исходный полутоновой снимок (256х256); б — гистограмма 


На рис. 3 представлены результаты 
кластеризации снимка, приведенного на 
рис.1 а базовым (рис. 3 а) и модифи- 
цированным (рис. 3 6) алгоритмами $ЕСМ. 
Несмотря на то, что в обоих случаях 
происходит выделение объектов интереса, 
предложенный метод позволяет обеспечить 
лучшее визуальное выделение области 
влияния гематомы за счет отсутствия 
излишней детализации фона. 
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На рис. 4 представлены результаты 
кластеризации снимка, приведенного на рис. 
2 а базовым (рис. 4 а) имодифицированным 
(рис. 4 6) алгоритмами $ЕСМ. Применение 
модифицированного алгоритма позволило 
получить качественное выделение различ- 
ных областей изображения при уменьшении 
излишней детализации по сравнению с 
применением базового метода. 
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Рис. 3. Кластеризация изображения рис. Та: а — базовым алгоритмом $ЕСМ, 
б — модифицированным алгоритмом $ЕСМ 


б 


Рис. 4. Кластеризация изображения рис. 2 а: а — базовым алгоритмом $ЕСМ, 
б — модифицированным алгоритмом $ЕСМ 


Выводы 

Предложенный в данной работе 
модифицированный алгоритм $ЕСМ позво- 
ляет обеспечить лучшее выделение объек- 
тов интереса и структуры изображения в 
целом по сравнению с базовым методом 
при выполнении кластеризации полуто- 
новых снимков. При этом сингулярное 
разложение — не единственный метод, 
который может быть применен при 
динамическом преобразовании исходных 
данных. Перспективным подходом, напри- 
мер, является использование для этой цели 
метода независимых компонент. Предло- 
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женный подход, предполагающий динами- 
ческое преобразование исходных данных, 
может быть применен также для других 


алгоритмов нечеткой — кластеризации, 
базирующихся на ЕСМ. 
Однако, следует отметить, что 


модифицированный метод обладает нес- 
колько меньшим быстродействием, т.к. при 
выполнении динамического преобразова- 
ния исходных данных требуется повторное 
выполнение кластеризации сетью Кохо- 
нена и пересчет матрицы расстояний и 
нечеткой функции принадлежности. Также 
модифицированный алгоритм в данной 
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реализации не может обрабатывать полуто- 
новые изображения в случае расширения 
пространства входных признаков. Для не- 
которых изображений условие (7) сраба- 
тывает на 1-ой итерации. В этом случае 
результаты кластеризации модифициро- 
ванного алгоритма ЕСМ будут такими же, 
как для базового метода. 
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Стаузсае ппасез Ё177у са$егтте 
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Дупапис ша! дайа тапзюгтайноп Базе оп 
зшеч[аг десотроз!оп оЁ ша! ааа апа Ра7ту 
тетбег$Ьр Рапсйоп (\ЫсВ 1$ пиегргае4 а$ 
ши[ааптеп$юпа! плазе) оп еасв чер оЁ 
сасшайоп. АЙег зшеи]аг десотроз1оп 110$ 
шЮгтайуе соати$ оЁ ей зшечаг уесюг$ 
таблх аге сВозеп  ааютайсаПу. Ты$ 
ашюотае соитиз ппрокапсе деесйоп 1$ 
Базе4 оп соеЁйслеп$ \ысВ аге саси]аед ув 
ие этешаг уесбюог$ тайлх изасе. То ауо1а 
пегайуе  шИчепсе оЁ зшещаг уаще 
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ргорозе арргоасб тау Бе гееазе4 Тог апу 
Ва7ху саегше аеогирт Базед оп ЕСМ. 

Тре ехрегитеша! гезеагсВе$ \уеге та4е 
юг уамои$ геа| отаузсае теФса| птазез. 
Зертешайоп гези $ оЁ зисВ итазез ш сотраге 
Уф Базе 3ЕСМ штефо4 аге зВо\уп ш 5 
агасе. ТБезе гезиз  Четопугае Фе 
еНесйуепез$ оЁ ргорозе4 сазегте ао т. 

Ргорозеа сачегие шефо Каз пех 
Ч1задуапасез: десгеазте регГюгтапсе ие ю 
гереайиех оГ Копопеп’$ шар сауегие Юг 
сетег$ ОР №77у сазег$ оМаште, Ч$апсе 
тайлх апа Ку  тшетфегзЫр  РЮпсйоп 
сасйайоп айег 4упаписа| фшапзЮгтайоп оЁ 
ша] аа; абзепсе оЁ ау ю \оЖ ми 
епрапсе аппепзюоп о шогтайуе аби щез Юг 
ргосеззе этаузсае ппазез; юг зоте птазез 
Дупаписа| ша! даа бапзЮгтайоп \Ш Бе 
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